Non-Markovian effect on the quantum discord 
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We study the non-Markovian effect on the dynamics of the quantum discord by exactly solving a 
model consisting of two independent qubits subject to two zero-temperature non-Markovian reser- 
voirs, respectively. Considering the two qubits initially prepared in Bell-like or extended Werner-like 
states, we show that there is no occurrence of the sudden death, but only instantaneous disappear- 
ance of the quantum discord at some time points, in comparison to the entanglement sudden death 
in the same range of the parameters of interest. It implies that the quantum discord is more useful 
than the entanglement to describe quantum correlation involved in quantum systems. 
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Since the proposal by E. Schrodinger entanglement 
has triggered off many famous and imaginative discus- 
sions for our deeper understanding of the quantum world. 
Nowadays, entanglement is considered to be not only a 
vital concept in physics but also a prime resource for 
quantum information processing (QIP) Q. 

Entanglement is something unique without classical 
counterpart, which is reflected in the implementation of 
quantum teleportation Q , quantum cryptography and 
universal quantum computing However, there seem 
some exceptions. For example, in Grover search 0, @] 
and in deterministic quantum computation with one pure 
qubit [8|, entanglement seems unnecessary in the imple- 
mentation, although the accomplishment of those quan- 
tum tasks results in something unaccessible in a purely 
classical way. This could be understood as that entan- 
glement only represents a special kind, but not all, of the 
quantum correlation in those systems. 

The present work focusses on another concept for 
quantum correlation, termed quantum discord [9l. [loL fTll| . 
which was introduced as the difference between two nat- 
ural quantum extensions of the classical mutual infor- 
mation. We have noticed some investigations about the 
quantum discord to work for quantum algorithms which 
showed that quantum discord is more practical than en- 
tanglement to describe quantum correlation fl2l . Il3j . Be- 
sides, quantum discord could be used to improve the ef- 
ficiency of the quantum Carnot engine and to better 
understand the quantum phase transition and the pro- 
cess of Grover search [l5|, [l(| . 

It is considered that the quantumness captured by dis- 
cord is different from entanglement [uj. This was 
also investigated under the Markovian environment in 
a recent publication for dissipative dynamics [l7| that 
the discord with an asymptotical decrease is more robust 
than the entanglement with sudden death under the same 
conditions. However, does this characteristic hold under 
non-Markovian environmental conditions ? 



'Electronic address: 
t Electronic address: 



bowangphysics@gmail.com 



mangf eng@ wipm .ac.cn 



Since the non-Markovian environment keeps memory, 
which could influence the dynamics of the coupled sys- 
tem, it would be of great interest in an investigation of 
the dissipative dynamics of the discord under the non- 
Markovian environment. To this end, we will make com- 
parison with 'entanglement sudden death' (ESD) [l~8l |. 
a widely studied terminology reflecting the fragility and 
complexity of the entanglement. There have been a lot 
of investigations on the ESD [H, H M, 
HE H3, HI] j including some experimental observations for 
ESD with the photon pairs and the atomic ensembles 
[2^ |. Meanwhile, the Markovian environment will also be 
treated here as a comparison. Our question is whether 
the quantum discord would present a similar behavior 
to the entanglement in the evolution with respect to the 
same characteristic parameters of interest. 

Our model consists of two independent qubits interact- 
ing, respectively, with their own reservoirs. To be simpli- 
fied, we assume the reservoirs to be in zero temperature. 
As the two qubits are initially entangled, we will show 
the dissipative dynamics of the quantum correlation re- 
flected in the time evolution of the quantum discord in 
comparison with the entanglement. 

We first present a brief review of the quantum dis- 
cord. In classical information theory, the Shannon en- 
tropy H(X) = -J2 x Px=x\°EPx=x is used to mea- 
sure the uncertainty of a random variable X, where 
Px=x is the probability with X being x. Similarly, 
the joint entropy, which measures the total uncer- 
tainty of a pair of random variables X and Y, is de- 
fined as H(X, Y) = -J2x,yPX=x,Y=yl0gPx=x,Y=y,wtiil 

Px=x,Y=y being the probability in the case of X — x and 
Y = y. As a result, the mutual information for the corre- 
lation between two random variables X and Y is defined 
as I(X : Y) = H(X) + H(Y)-H(X,Y), whose quantum 
version can be written as, 

T(X :Y) = S(px) + S(py)-S(pxy), U) 

where S(p) = — Tr^logp) is the von Neumann entropy 
of p, and px(py) is the reduced density matrix of pxY 
by tracing out Y(X). 

For classical probability distributions, the Bayes rule 
Px\Y=y = Px,Y=y/PY=y leads to an equivalent expression 
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for the mutual information 



I(X : Y) = H(X) - H(X\Y), 



(2) 



where H(X\Y) = J2 y PY= y H(X\Y = y) = 

- Ys X , y Px=x,Y=y \ogPx=x\Y= y is the conditional entropy 
of the random variables X and Y for the average uncer- 
tainty about the value of X given that the value of Y 
is known. In order to generalize Eq. (2), we measure 
the subsystem Y by a complete set of projectors {Hi}, 
corresponding to the outcome i, which yields px\i = 
TrY(ILipxY^h)/pi, with pi = TrxY^iPxY^-i)- 
So we may define the quantum conditional entropy 
S {ni} (X\Y) = T,iPiS(px\i)- Following Eq. ©, we have 
the quantum mutual information alternatively defined by 



J {Ui} {X :Y) = S(p x ) - S { n t} (X\Y). 



(3) 



The above quantity strongly depends on the choice of the 
measurements {IE}. By maximizing i7m < }(-X' : Y) over 
all {IE}, we define an independent quantity J{X : Y) = 
max {ni }{S{px) - EiPiS(px\i)} = S(p x ) - S(X\Y) as 
a measure of the classical correlation. 

Having the two quantum analogs of the classical mu- 
tual information T(X : Y) and J(X : Y), we define the 
difference between them as the quantum discord, 



V(X : Y) = T(X : Y) - J(X : Y), 



(4) 



reservoir [291 ]. In our model, we use the Lorentzian 
spectral distribution 



J(cj) = 



1 



7oA 2 



2tt {lo - lu) 2 + A 2 



(7) 



where the parameter A, defining the spectral width of the 
coupling, is connected to the reservoir correlation time 
tb by the relation tb ~ A -1 . For our purpose, we define 
another parameter 70 regarding the decay of the atomic 
excitation in the Markovian limit of the flat spectrum. 
The relaxation time scale tr over which the state of the 
system changes is then related to 70 by tr w 7^ . The 
Markovian and the non-Markovian regimes are distin- 
guished by the relation of the parameters 70 and A. In the 
Markovian regime there is 70 < A/2 or tr > 2tb, and the 
non-Markovian regime corresponds to 70 > A/2 or tr < 
2tb [13, H^. In the non-Markovian regime, straightfor- 
ward solution of Eq. © could yield d{t) = Ci(0)x(t) 
where x(*) = e~ xt / 2 [cos {dt/2) + (X/d) sin (dt/2)] with 
d = ^/|27oA — A 2 |. In contrast, in the Markovian regime, 
C\{t) has the similar form but with cos[-](sin[-]) replaced 
by cosh[-](sinh[-]). So the dynamics of a qubit S can be 
represented by the reduced density matrix 



p S it) 



pfi(o)xW 2 pfoWxW 

Poi(0)x(*) l-pfi(0)x(t) 2 



(8) 



which is interpreted as a measure of the quantum corre- 
lation [ItMllll- 

For two noninteracting qubits, i.e., A and B (involving 
two levels in each) locally interacting with the reservoirs 
Ra and Rb, respectively, the Hamiltonian of each sub- 
system (i.e., qubit+reservoir) is given by 



H 



LJq(T + <J- 



-^2(g k b k a++g k \bla-), (5) 

k 



where loq denotes the transition frequency of the two- 
level system (i.e., the qubit) with <r± the correspond- 
ing atomic raising and lowering operators. The index 
k labels different field modes of the reservoir with fre- 
quencies uj k . bf. (6fe) is the creation (annihilation) oper- 
ator of the reservoir field with g k the coupling constant 
to the qubit [2(| H3, H^]. We assume that the initial 
state of the qubit with the zero-temperature reservoir 
is I* {0)) § = (C o (0)|0)s + Ci(0)|l) s ) |0) fls . The ampli- 
tudes at any time can be obtained exactly by, 



Ci(t) 



C (t) = C (0), 
dtFty-t^dtf), 



(6) 



where the correlation function F(t — t ) = 

/(LJ(w)e , ( w "- ,J ) ( '- f ' with J(uS) the spectral den- 
sity of the reservoir. The exact form of C\(t) depends 
on the particular choice of the spectral density of the 



where pfj(0) = Ci(0)C*(0), x( f ) 2 has discrete zeros at 
t n = 2 [nn — arctan(d/A)] /d, with n being an arbitrary 
integer. Because the system is composed of two noninter- 
acting parts evolving independently according to Eq. |5]), 
we may simply construct the density matrix p for the 
two-qubit system by the reduced single-qubit density ma- 
trices [H Hi). In the basis {|11), |10), |01), |00)}, we 
measure the qubit B from the matrix p(t) by projecting 
on {cos6>|l) B + e^sin0|O) B ,e-^sin0|l) B - cos(9|0) B }. 
Then the quantum discord could be calculated numeri- 
cally using Eq. Q . 

The initial states we assume are the extended Werner- 
like states (EWL) [13, [S3], defined as 



(9) 



where |f) = |$) or |*) with |$) = a|10) + (1 - a 2 ) 1 / 2 ^)!) 
and I*) = a 1 00) + (1 - a 2 )^ 2 \ll). For r = 0, the EWL 
states become totally mixed, while they reduce to the 
Bell-like pure state |$) or |W) in the case of r = 1. 

In what follows, we will check if sudden death happens 
in the dynamics of the discord under the same condition 
with respect to the entanglement. To this end, we em- 



ploy concurrence representing the entanglement 31| as a 
comparison. We get started from the Markovian regime. 
In Fig. 1, we plot the dynamics of the concurrence and 
the discord for the same Bell-like state |\&) as a function 
of the dimensionless quantities j^t and a 2 , under the 
same condition A/70 = 10. There are clearly two ranges 
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FIG. 1: (Color online) Variation of the concurrence (a) and 
the discord (b) from the initial state |$) with respect to the 
dimensionless quantities 70 i and a 2 , in the Markovian case 
(e.g., A/70 = 10). 





FIG. 2: (Color online) Variation of the concurrence (a) and 
the discord (b) from the initial state |$) with respect to the 
dimensionless quantities jot and a 2 in the non-Markovian case 
(e.g., A/70 = 0.1). 



for the variance of the concurrence: For a 2 < 1/2, the 
ESD occurs, which denotes the entanglement vanishing 
abruptly after a finite time; But for a > 1/2, the concur- 
rence vanishes in an asymptotical way. In contrast, the 
discord vanishes only asymptotically with the variance of 
those characteristic parameters. 

Fig. 2 presents the case for the non-Markovian regime 
(i.e., A/70 = 0.1) from the initial Bell-like state |'F). It 
can be seen that the discord of |^) periodically vanishes 
in accordance with the zero points of the function x(t) 
following the asymptotical damping. On the contrary, 
there is ESD for the concurrence under the same con- 
dition in some ranges of the parameter a 2 . A clearer 
difference between the behaviors of the discord and the 
concurrence could be found in Fig. 3(a). 

To see what happens in the mixed state, we have made 
a study in Fig. 4 for two initial EWL states p% WL and 
Pewl w ith a — 1/ V^2 in the non-Markovian regime. We 
have found that the concurrence exists at the beginning 
of the evolution (i.e., t = 0) only for r > 1/3, while 
the discord is always positive for r > 0. The common 
feature of the dynamics for both the concurrence and the 
discord is the decrease with the decrease of the purity r, 
implying that the mixedness of the initial state affects 
both of them. The difference between them is similar to 
the case of the pure state, as shown more clearly in Fig. 
3(b). 

How much the non-Markovian environment influences 
the discord? Since A/70 characterizes the degree of the 
non-Markovian effects, we have studied the dynamics of 
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FIG. 3: (Color online) The comparison of the evolution of 
the concurrence with the discord in the non-Markovian case 
(e.g., A/70 = 0.1): (a) For the initial state |^) and (b) for 
the initial mixed state p% WL . The green dashed line denotes 
the discord. The black (dotted and solid) curves correspond, 
respectively, to the concurrence with different parameters a 2 
or r. o denotes the zero points of the discord. 





FIG. 4: (Color online) Variation in the non-Markovian case 
with respect to the dimensionless quantities 7oi and r, where 
A/70 = 0.1 and a 2 — 1/2, (a) and (b) are for the concurrence 
and the discord, respectively, from the initial states p%wl\ 
(c) and (d) correspond to the same variables as in (a) and 
(b), but from the initial states p* wi- 



the discord in Fig. 5 starting from the EWL states Pewl 
with r = 1 and a = l/v3. We can observe that the 
revival amplitude of the discord increases with the de- 
crease of the A/70 corresponding to the enhancement of 
the non-Markovian effects. 

Our results for the non-Markovian regime have shown 
something very different from the Markovian regime [l7| : 
The discord disappears instantaneously at some time 
points following the asymptotical dissipation. In some 
sense, this could also be called sudden death and sudden 
revival of the discord, but we prefer to call it the instan- 
taneous disappearance. Although the mechanism behind 
the instantaneous disappearance needs further clarifica- 
tion, some points seem clear to us: The revival is due to 
the memory effect of the non-Markovian reservoir. The 
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FIG. 5: (Color online) Variation of the concurrence (a) and 
the discord (b) from the initial states p%wl with respect 
to the dimensionless quantities 70 1 and A/70 in the non- 
Markovian case, where the parameters r = 1 and a 2 = 1/3 
are set in the calculation. 



change of A in Fig. 5 implies that the stronger memory 
effect of the reservoir yields the larger revival amplitude 
of the discord. Besides, the zero discord in our model 
happens in the case of the evolution to the component 
state 1 00), 1 1} , 1 10) , |11) or the maximal mixed state 
I. In particular, for some values of the parameters, such 
as a = or 1 in Figs. 1 and 2, and r — in Fig. 4, 
the discord remains zero in the evolution. In contrast, 
for other case, the discord only falls asymptotically, with 
some fluctuations reaching zero instantaneously. 

The Lorentzian spectral distribution in Eq . (7) had 
been widely employed in quantum optics [291 ] . and was 
recently used in QIP studies [20]. Since strong coupling 
between matter and li ght has been available experimen- 
tally in some systems [32|], dynamics on short-time scale 
could in principle be observed and even be manipulated. 
As a result, non-Markovian effect would become more 
and more important in the exploration of QIP under 
real experimental environment. Moreover, our studies 



for both Markovian and non-Markovian reservoirs have 
clearly demonstrated that the discord is more robust in 
dissipative evolution with respect to entanglement and 
more suitable to describe the quantum correlation in- 
volved in the system. As a result, it is understandable 
why some quantum algorithms could work well in the ab- 
sence of entanglement. The discord should be larger than 
zero in the implementation of those algorithms in the case 
that the entanglement approaches zero. A deeper study 
in this respect is expected. 

In conclusion, we have investigated the dynamics of the 
discord using an exactly solvable model where each qubit 
independently interacts with its own zero-temperature 
reservoir. We have compared the discord with the con- 
currence using the same initial states and the same reser- 
voir conditions. We have also discussed the different ef- 
fects from the Markovian and non-Markovian reservoirs. 
Further work would be to explore the dynamics of the 
discord subject to some specific decoherence mechanism, 
and to quantitatively evaluate the necessary discord re- 
quired in implementing different quantum algorithms. 
We will also consider the situation at finite temperature. 
In comparison with some recent work on non-Markovian 
effect on the dissipation of the system in a microscopic 
way [HI, our present work from the phenomenological 
viewpoint, might be more practical to explain some ex- 
perimental observations of the dissipation of quantum 
correlation subject to a realistic environment. 

This work is supported by NNSFC under No. 
10774163, No. 10774042 and No. 10775175. 

Note added: The work was submitted to Physical Re- 
view A in September. But we have just been aware of a 
similar work in arXiv:0911.1096v2. 
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